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Abstract 

The two-dimensional steady-state Boltzmann equation for hard-disk molecules in the presence of 
a temperature gradient has been solved explicitly to second order in density and the temperature 
gradient. The two-dimensional equation of state and some physical quantities are calculated from 
it and compared with those for the two-dimensional steady-state Bhatnagar-Gross-Krook(BGK) 
equation and information theory. We have found that the same kind of qualitative differences as 
the three-dimensional case among these theories still appear in the two-dimensional case. 
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I. INTRODUCTION 



The behaviors of gases in nonequihbrium states have received considerable attention from 



the standpoint of understanding the characteristics of nonequihbrium phenomena. 
1^,0,0,0, ^ The kinetic theory has contributed not only to the understanding of nonequihb- 
rium transport phenomena in gases but also to the development of general nonequihbrium 
statistical physics. It is well accepted that the Boltzmann equation is one of the most re- 
liable kinetic models for describing nonequihbrium phenomena in gas phases. In the early 
stage of studies on the kinetic theory, great effort has been paid for solving kinetic model 
equations such as the Boltzmann equation and deriving nonequihbrium velocity distribu- 
tion functions and macroscopic nonequihbrium transport equations in terms of microscopic 
molecular quantities. These attempts were strongly related to the development of general 
nonequihbrium statistical physics such as linear no neq uihbrium thermodynamics, Onsager's 
reciprocal theorem and the linear response theory. [lO|, |ll | 

Among various methods which give normal solutions of the Boltzmann equation, the 
Chapman-Enskog method has been widely accepted as the most reliable method. It had 
been believed that Burnett determined the complete second-order solution of the Boltz- 
mann equation by the Chapman-Enskog method. Physical importance of the 
second-order coefficients has been also demonstrated for descriptions of shock wave profiles 
and sound propagation phenomena, jl^. Iisl [3| However, it was reported that Burnett's 
solution is not complete, and Schamberg derived the precise velocitv distribution function 
of the Boltzmann equation to second order for Maxwell molecules. [l7| On the other hand, 
because of its mathematical difficulty, the complete second-order solution of the Boltzmann 
equation for hard-core molecules has been derived quite recently. Its validity has been 



also demonstrated by numerical experiments of both a molecular dynamics simulation and 



a direct simulation monte carlo method. 
Gross-Krook (BGK) equation I21I 



22, Ul 



19L 
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Other kinetic models like the Bhatnagar- 

H H 22, Q have been proposed mainly 



to avoid the mathematical difficulties in dealing with the collision term of the Boltzmann 
equation. Although it was believed that its results approximately agreed with those of the 
Boltzmann equation, some qualitative disagreements have been found in the second-order 
solutions. |3] 

Following its success and usefulness, the Boltzmann equation is widely used in order 



2 



30 



31 



to describe vario us g as-phase transport _2henomena such as granular gases j2 
plasma gases 0, Is^ . polyatomic gases 



35l |. relativistic gases |36l| and chemically reacting 
gases. Chemical reactions in gas phases have been studied with the aid of gas collision theory. 



For the differential cross-section of chemical reaction, the line-of-centers model proposed 
Present has been accepted as a standard model to describe the chemical reaction in gases. 
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by 



4J| This model can be derived explicitly using a collision law of hard- 



core molecules, and the diameter of hard-core molecules is regarded as a distance between 



centers of monatomic molecules at contact. 



371 . l38l l42| In equihbrium states, experimental 



results includirig the t emp erature dependency can be fitted by the results from the line-of- 
centers model. I37I. [4^ I41I Under nonequilibrium situations such as gases under a heat flux 
or a shear flow, their pure nonequilibrium contributions to the rate of chemical reaction 

39, 43, 44, 45, Q, 47, Q, Q 



have also attracted much attention. 



Since nonequilibrium 



correction terms of the chemical reaction rate are quadratic functions of nonequilibrium 
fluxes, the explicit nonequilibrium velocity distribution function of the Boltzmann equation 
for hard-core molecules to second order is needed to derive it based on the line-of-centers 



.0,0,0,0 



model. The pure effect of a heat flux on the chemical reaction rate has been 

recently calculated using the second-order velocity distribution function of the Boltzmann 
equation for hard-core molecules. In the letter, a thermometer to measure a relation 
between a kinetic temperature of gases under a heat flux and a temperature of a heat bath 
has been also proposed. 

It is one of the most significant subjects in modern statistical physics to construct the 
nonlinear nonequilibrium statistical mechanics and thermodynamics for a strongly nonequi- 
librium state beyond the local equilibrium state, called the local nonequilibrium state. 
Zubarev0. Isi] has developed nonequilibrium statistical mechanics and obtained the gen- 
eral form of a nonequilibrium velocity distribution function with the aid of the maximum 
entropy principle. It is expanded to first order under some constraints to obtain the first- 
order nonequilibrium velocity distribution function. 0| Jou and his coworkers have derived 
the nonequilibrium velocity distribution function to second order by expanding it to second 



53, 



54| Information theory 



order under some constraints, which is called information theory, 
has attracted interest in the development of a general framework for nonlinear nonequilib- 
rium statistical mechanics which can describe the local nonequilibrium state. The nonequi- 
librium velocity distribution function from information theory has been applied to nonequi- 
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librium systems, and some predictions were made. For example, in dilute gas systems under 
nonequilibrium fluxes, an anisotropic pressure and a nonequilibrium ternperature which is 
not identical with the kinetic temperature have been predicted. 55|, |56|, |57|, 15811 There are 
also several applications of information theory to chemically reacting gases. [3^, |4^ How- 
ever, some qualitative differences between information theory and the Boltzmann equation 
have been recently reported, and the invalidity of information theory as universal nonlinear 
nonequilibrium statistical mechanics has been demonstrated, jsof 

There have been no efforts to solve two-dimensional kinetic models to second order 
and discuss the two-dimensional second-order nonequilibrium phenomena, thou gh two 
dimensional transport phenomena have created great interest. j60|, |6l|, |62, |63|, l64L 
The main aims of this paper are to reconstruct all the results obtained in refs. |l8l. 
in the case of two dimensional, and to discuss properties of the two-dimensional nonlinear 
nonequilibrium phenomena which reflect the local nonequilibrium state. In Secs|n] and IIIH 
we have derived the explicit velocity distribution function of the two-dimensional steady- 
state Boltzmann equation for hard-disk molecules to second order by the Chapman-Enskog 
method. In order to achieve that, we have extended the method we developed in ref. jl^ to 
the two-dimensional case. We also obtain the nonequilibrium velocity distribution functions 
to second order for the two-dimensional steady-state BGK equation and information theory 
in Secs ilV and IIVBI respectively. All the nonlinear nonequilibrium velocity distribution 
functions are graphically compared in SecEl Using the two-dimensional nonequilibrium 
velocity distribution functions to second order, we discuss differences among those theories 
appearing in the two-dimensional nonlinear nonequilibrium transport phenomena in Sec I VII 
In Sec iVIH we explain how to calculate the effect of steady heat flux on the rate of chemical 
reaction based on the line-of-centers model in the two-dimensional case, and apply the two- 
dimensional nonequilibrium velocity distribution functions to second order to calculate it. 
We have also investigated dimensional dependency appearing in the nonlinear nonequilib- 
rium phenomena which reflect the local nonequilibrium state. Our discussion and conclusion 
are given in Sec JVIIll 
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II. THE CHAPMAN-ENSKOG METHOD FOR SOLVING THE TWO- 
DIMENSIONAL STEADY-STATE BOLTZMANN EQUATION 



Let us introduce the Chapman-Enskog method to solve the two-dimensional steady-state 
Boltzmann equation in this section. Assume that we have a system of dilute gases in a 
steady state, with velocity distribution function /i = /(r, vi). The appropriate steady-state 
Boltzmann equation is 

Vl-V/i = J(/i,/2), (1) 

where the collision integral /2) is expressed as 

J(A, f,)^l I - /i/2)^/dWv2, (2) 

with f[ = /(r, v']^) and = /(r, Vg): v'^^ and are postcollisional velocities of Vi and 
V2, respectively. The relative velocity of two molecules before and after an interaction 
has the same magnitude g = |vi — V2I ; the angle between the directions of the relative 
velocity before and after the interaction is represented by x- The relative position of the 
two molecules is represented by b, called the impact parameter, (see Fig^ The impact 
parameter b depends on kinds of the interaction between molecules, and one should specify 
the intermolecular interaction so as to explicitly determine the impact parameter b in the 
collision term J(/i, /2). Note that x can be expressed as a function of b for a central force. 
Suppose that the velocity distribution function fi can be expanded as 

/i = fi'' + fi'' + ff + ■■■ = fi'\l + + 0f ^ + ■■■), (3) 

where the small expansion parameter will turn out the Knudsen number K = l/L, which 
means that the mean free path of molecules I should be much less than the characteristic 
length L for changes in macroscopic variables, fi^^ is the local Maxwellian distribution 



function, written as 



(o)_ n[r)m 



2/tT(r) 

with m mass of the molecules and k the Boltzmann constant. n{Y) and T(r) will be identified 
later as the density and the temperature at position r, respectively. Substituting eq.© into 
the two-dimensional steady-state Boltzmann equation (Q), we arrive at the following set of 
equations which we will solve completely in this paper: 

= VI ■ V/f \ (5) 



to first order and 

^[/n#=vi-v/«-j(/«,/«), (6) 

to second order. The linear integral operator is defined as 

)]Xi = /f ^ {X[ -X,+X',- X,)gAhdw,. (7) 

The solubility conditions of the integral equation (0) are given by 

j $iVi • V/J°)dvi = 0, (8) 

where $i are coUisional invariants: 

$1 = 1, $2 = mvi, $3 = ^mv^. (9) 

Substituting eq.Q into the solubility conditions (jH)), it is seen that n/tT is uniform in the 
steady state. We use this result in our calculation to second order. Similarly, the solubility 
conditions of the integral equation (jH)) are given by 

$iVi ■ V/i^'^dvi = 0, (10) 

which will be considered in Sec illl C"2l 

To construct solutions of the integral equations (0) and © definite, four further conditions 
must be specified; we identify the density: 

n(r) = j /idvi = j /^dvi, (11) 

the temperature: 

n(r).r(r) ^ / !^Adv, = / ^/!»'dv„ (12) 

and the mean flow: 

Co = y mvi/idvi = j mvi/f ^dvi. (13) 

Here we assume that no mean flow, i.e. Co = 0, exists in the system. The introduction of 
these conditions distinguishes the Chapman-Enskog adopted here from the Hilbert method 
in which the conserved quantities are also expanded. 0] We assert that conditions (fTT|) . (fT^ 
and (fT^ do not affect all our results in this paper. It should be noted that, to solve the 
integral equations (0) and (jH)), we should consider only the case in which the right-hand 
sides of eqs.(jni) and (jHl) are not zero: if the right-hand sides of eqs.© and (jHl) are zero, 
the integral equations (0) and (0) become homogeneous equations which do not have any 
particular solutions. P] 



III. A METHOD FOR SOLVING THE INTEGRAL EQUATIONS 



A. A general form of the velocity distribution function 

To solve the integral equations (0) and We assume a general form of the velocity 
distribution function: 

fe oo 



/l = /i 



(0) 



r=0 k=l r=0 



(14) 



with Ci = (m/2KT)^/^vi the scaled velocity. S^{X) is a Sonine polynomial, defined by 

oo 

(1 - uj)-'-'e-^ = Y^r{p + k + l)Sl{X)u^, (15) 

p=0 

and 

Yfcr(Cl) = BkrYk{Ci) + CkrZkiCi), (16) 

where B^r and Ckr are coefficients to be determined. Introducing the polar coordinate 
representation for Ci, i.e. c\x = ci cos 6*, ciy = ci sin^, 

Yk{ci) = 

\ m J 

and 



n(ci) = cfcosfc0, (17) 

\ m / 



k 



c\smk(t). (18) 

Assumption of the velocity distribution function of eq. (jl4p has some mathematical ad- 
vantages in our calculation. Firstly, it is sufficient to determine the coefficients Bkr-, because 
Ckr can be always determined from Bkr by transformations of axes. Secondly, some impor- 
tant physical quantities are related to coefficients Bkr and Ckr'- e-d- the density (fTTj) . the 
temperature (fT^ and the zero mean flow (fT^ with /i in eq. (fT^ lead to the four equivalent 
conditions: 

5oo = l, Eio = Cio = 0, 5oi = 0. (19) 
Similarly, the pressure tensor Pij deflned by 

Pij = I ) / dcimciiCij/i, (20) 



for i,j = X and y is related to B20 and C2o- 

The coefficients Bkr except for those in eq. ()19|) can be calculated as follows. Multiplying 
the two-dimensional steady-state Boltzmann equation (H)) by 



Qfer(ci, - . ^2«:T; r(A; + r + l) 
and then integrating over (2KT/m)^/^Ci, it is found that 



\ m 



< Cl ■ VQkr > 



/2kT 
\ m 



c 



2kT 



m 



I I j^Q'kr-Qkr)hf29dhdC2dc^, (22) 



where r(X) is the Gamma function, < X >av= {^kT Im) J X/idci and represents the 
postcoUisional Qkr- We should calculate both sides of eq.(j22I) for every k and r, because 
eq. (j22p leads to simultaneous equations to determine B^r- 

For convenience, we introduce Qkr and A^r as the left-hand and the right-hand sides of 
eq.(I221), respectively: 



Qkr — 



2kT 



m 



< Cl ■ VQkr >av +V 



\ m 



< CiQkr > c 



(23) 



and 



A 



kr 



2kT 



m 



III ~ Qfcr)/l/2^dWC2dCi. 



(24) 



We will calculate Qkr and Aj^,. separately. From Appendix ^ the result of Qkr becomes 

'2kT\^ 



m 



. , k 1 d.,T dyT, 

('^ + 2 - 2^iDk,r^ + Ek,r^] 



k,r—l' 



+ E, 



k,r—l 



d T 



'x-l^k,r 



dyEk. 



(25) 



where diX denotes dX/di for i = a; and y. Dk^r and Ek^r are functions of Bkr and Ckr, as is 
written in Appendix 1X1 



B. The collision term A 



kr 



Next we calculate the collision term A^^ in eq. ()24|) . We should specify the kind of the 
interaction of molecules so as to perform the calculation of the collision term A^^- For 
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hard-disk molecules, the impact parameter b is given by the relation 

6 = dcos|, (26) 

where d is the hard-disk molecular diameter. The collision differential cross section is ob- 
tained by 

d6 = -^sin|rfx, (27) 
Therefore, A^r for hard-disk molecules, A^^, becomes 

A^. = ^ I'^FLix) - sin |dx, (28) 

where Fj^^{x) is defined as 

Kix) ^ / / Q'krflf29''dC2dc„ (29) 

and we have used ^^^^(0) = {2KT/mf J J Qkr f 1/29^ dc2dci. Note that x = if 6 > 

From eq. ipHj) . it is sufficient to calculate F^^{x) ^oi deriving A^^. The details of F^^^x) ci-re 
written in Appendix IB II Several explicit forms of A^^ are also demonstrated in Appendices 
O and El From the definitions ()23|1 and (j2l|l . both sides of eq. ()22j) for arbitrary k and r can 
be calculated for hard-disk molecules via 

= AH , (30) 

which produces a set of simultaneous equations determining the coefficients Bkr, as is ex- 
plained in Appendices O and O Here denotes Qkr for hard-disk molecules. 



C. Determination of Bkr 

We will determine the first-order coefficients Bl^. and the second-order coefficients B^^ 
in accordance with the previous two subsections, which corresponds to solving the integral 
equations (0) and ©, respectively. Here the upper suffices I and II are introduced to specify 
the order of K. 
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1. The First Order 



We show the results of the first-order coefficients of which the solution of the integral 
equation \ is composed. They can be written in the form: 

2d T 

Bkr = 5k,ibir—7=——;- (31) 
y iTidni 



the Knudsen number K. Though B],^ was 
6^ . i.e. B\^ for r = 1, we have obtained 



Values of the constants hir are given in Table HI The calculation of B\^ is explained in 
Appendix O It is seen that B\^ is of the order o: 
derived only to the lowest order approximation 
B\^ for r < 7 in this paper. This ensures the convergence of all the physical quantities 
which will be calculated in this paper. It should be mentioned that our value of B} for the 
lowest Sonine approximation, i.e. r = 1, is identical with Sengers's value |62j|. Once Bl^ have 
been calculated, Cl^ can be written down directly by replacing dxT by dyT by symmetry. 
Substituting all the first-order coefficients derived here into eq. ljlfjj) . we can finally obtain 
the first-order velocity distribution function f[^\ 

(2) 

2. Solubility Conditions for (p^ 

Since the first-order velocity distribution function f[^^ has been obtained, the solubility 
conditions for the integral equation (jH)) should be considered before we attempt to derive 
the explicit second-order solution The solubility conditions for (f)f\ eqs. ^TU^ . lead to 
the condition 

V • J(^^ = 0, (32) 
where J*-^-*, i.e. the heat fiux for fi^\ can be obtained as 



._,„(^)'^Vr. (34) 

with the appropriate value for hu listed in Table |l] It must be emphasized that, since J^^-*, 
i.e. the heat fiux for fl , does not appear, the solubility conditions of the two-dimensional 

(2) 

steady-state Boltzmann equation for 0^ lead to the heat fiux being constant to second 
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order. This fact is in harmony with a general property that the total heat flux should be 
uniform in the steady state. From eqs.()32|) and ()34|). we also obtain an important relation 
between {VTf and V^T, 

^ + V^T = 0. (35) 
Owing to the relation ()35|) . terms of V^T can be replaced by terms of (VT)^. 

3. The Second Order 

We write down the results of the second-order coefficients B^f}^. of which (f)f^ is composed. 
Using the relation (j35|l . we can determine the second-order coefficients -BjJ. appearing in 
eq.dH)) as 

(36) 

Values for the constants 6or are summarized in Table IHl The calculation of Bj^^ is shown in 
Appendix ini We have calculated to 7th approximation, i.e. B^^ for r < 6 in this paper. 
The other second-order coefficients Bj}^ in eq. ()l(j|l can be written in the final form: 

Values for the constants b^^ and 6^,, are summarized in Table IIIIl The calculation of B]^^ 
is explained in Appendix |D| Owing to the properties of the trigonometrical function, C].|. 
can be obtained by replacing {d^Tf - {dyTf and d^T - d'^T in eq.§7^ by 2d^TdyT and 
2TdxdyT, respectively, using an axis change. 

One can see that both of BH and Bg are of the order of K^. As is also found in three 
dimension case|3], we have confirmed the fact that Bj}^ for = 4 and 6 do not appear. 
This fact stron gly suggests that B]^^ for all k greater than 2 do not appear, which was also 
expected in ref. jisl] and recently confirmed in ref.0|. We finally obtain f[^'^ by substituting 
the second-order coefficients obtained here into eqs. lfT^ and (fT^. Finally, we note that, 
though we have derived all the constants bir, bor, &2r in forms of fractions, we have 

written them in forms of four significant figures in this paper, since the forms of those 
fractions are too complicated. 
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4- The Velocity Distribution Function to Second Order 



The velocity distribution function for hard-disk molecules which we have derived in this 
subsection valid to second order is now applied to a nonequilibrium steady-state system 
under the temperature gradient along x-axis. In this case, the form of BqJ. in eq.pH|) becomes 

and, using the relation (j35p . in eq.(|37p can be transformed into a more simple form: 

(39) 

where values for the constants are summarized in Table HV1 The other second-order term 
becomes zero. 

From eqs. p4|) and ()16p . the velocity distribution function of the steady-state Boltzmann 
equation for hard-disk molecules to second order in the temperature gradient along x-axis 
can be written as 



biTUKT \2kT/ 



r>2 r>0 



(40) 



where the specific values for bir, feor and b2r are found in Tables HI ITTl and ll VI respectively. 
Jx corresponds to the x component of the heat flux in eq. ljMjl . Note that we have changed 
Ci to c. As can be seen from eq.()4U|). the explicit form of the velocity distribution function 
for hard-disk molecules becomes the sum of an infinite series of Sonine polynomials. 

Figure 121 gives the 0"^^) in eq. ipUl) scaled by mJ^/ri^K^T^ with the 3th, 4th, 5th, 6th and 
7th approximation 6or and 62r- It should be mentioned that, as FigO shows, the scaled cp^"^^ 
in eq. ()4()j] seems to converge to 7th approximation. Figure El provides the explicit form of 
the scaled 0*^^^ for hard-disk molecules with 7th approximation 6or and 62r- It is seen that 
the scaled 0^^-' for hard-disk molecules is strained symmetrically. 
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IV. OTHER NONEQUILIBRIUM VELOCITY DISTRIBUTION FUNCTIONS TO 
SECOND ORDER 



A. The Chapman-Enskog Solution of the two-dimensional Steady-State BGK 
Equation to second order 

For comparison, we also derive the velocity distribution function for the two-dimensional 



25 



26 



27 



steady-state BGK equation to second order by the Chapman-Enskog method. 
Suppose a nonequilibrium system subject to a temperature gradient along the x-axis in a 
steady state whose velocity distribution function is expressed as / = /(x, v). The steady- 
state BGK equation is written as 

v.dj = ll^E^^ (41) 
r 

where the relaxation time r dependent on position x through the density n{x) and 
the temperature T{x). Jle is the usual local equilibrium velocity distribution function 
fLE{x,v) = (mn(x)/27r/tT(x)) exp[— mv^/2/tT(x)]. It should be mentioned that, for the 
conservation laws, the collision term for the steady-state BGK equation, the right-hand side 
of eq.dH}, must satisfy 

j ^-jLEd^ = j <^ifd^r, (42) 

with the four collision invariants $i introduced in eq.Q. The velocity distribution function 
/ can be expanded as 

/ = /(°) + /(^) + f^'^ + ■ ■ ■ , (43) 

with /(°) = /le- Substituting eq. (jIH|) into the steady-state BGK equation (jSj), we arrive at 
the following set of equations: 

-^ = t;.9./W, (44) 

T 

to first order and 

-^ = v,dj^'\ (45) 

T 

to second order. It is found that equation ()44|1 with the requirement ()42|1 leads to nuT 
being uniform, and that equation (03)) with the requirement ()42|) leads to heat flux Jr^ 
calculated from eq.(pH|) being uniform. Using these facts, from egs. ^lj) and the velocity 

13 



distribution function to second order for the the two-dimensional steady-state BGK equation 
becomes 

n=0 

with the uniform heat flux = —2nH?TTdxT/m. 



B. Two-Dimensional Information Theory 



Let us construct two-dimensional information theory. [53|, |5J] The Zubarev form for the 
nonequilibrium velocity distribution function under a heat flux can be obtained by maxi- 
mizing the nonequilibrium entropy, defined as 



S{x)^-t, I /log /civ, (47) 
under the constraints of the density: 

n{x) = I fdy,, (48) 



and the temperature: 

n{x)KT{x) = J ^^fdv. (49) 
We assume no mean flow: 

j mvfdv = 0, (50) 
where denotes the zero vector. Furthermore, we adopt the heat flux as a constraint: 



= j -i^^xfdv. (51) 

It should be emphasized that nnT and the heat flux are now assumed to be uniform by 
contrast with the case for the steady-state Boltzmann equation where its solubility conditions 
lead to ukT and being constant to second order. 

We have finally derived the nonequilibrium velocity distribution function to second order 
in the heat flux by expanding the Zubarev's nonequilibrium velocity distribution function 
to second order as 

(52) 

14 



In eq. fl52j) . we have expanded the nonequilihrium temperature which has been obtained as 
O = T(l — SmJ^/An^K^T^) for two dimension. Such the modified velocity distribution 
function has been also obtained and used in three dimension. 13 note that, in all the 
macroscopic quantities calculated in this paper, there are no differences between the results 
from the modified velocity distribution function and Jou's velocity distribution function 
where nonequilihrium temperature is not expanded [ssl. Is^. Actually, the identifications of 
the density, the temperature and the mean flow in eqs. pi|) . ()12|) and ()13|) do not affect the 
physical properties of the velocity distribution function for the two-dimensional steady-state 
Boltzmann equation^ |, and those identifications must be satisfied for the conservation laws 
in the case for the two-dimensional steady-state BGK equation, (see eq.(jl2I)) 



V. DIRECT COMPARISON OF THE SCALED 

Figure E] exhibits the direct comparison of the scaled (p^'^^s for hard-disk molecules (j4m) to 
7th approximation with those for the steady-state BGK equation pUj) and information theory 
dSD). We have found that, as FigEl explicitly shows, the second-order velocity distribution 
function for hard-disk molecules ()40|) definitely differs from the others. We emphasize that 
such a difference never appears to first order. 



VI. NONLINEAR NONEQUILIBRIUM TRANSPORT PHENOMENA 

We can introduce the general form of the heat flux as 

J, = -wT^d^T, (53) 

where if indicates temperature dependence of the thermal conductivity and w \s a. constant 
that depends upon microscopic models. For example, if is calculated as 1/2 for hard-disk 
molecules, and w is determined as 2hiiK / d{K / -KmY^"^ with hn ~ 1.030 for hard-disk molecules 
(see eq. ()34|) ). Note that if and w cannot be determined explicitly from the BGK equation 
and information theory. From eq. ()53|) . the temperature profile T{x) in the nonequilihrium 
steady state can be determined as 

T(x) = [T(0)^+i-((^ + l)— x]^^. (54) 

w 

It is seen that the temperature profile T[x) becomes nonlinear except for 9? = 0. 
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Using eq. (pn|) . the equation of state in the nonequihbrium steady state can be obtained 

as 

777 

P., = nnT[5., + (55) 

with the unit tensor 5ij and the tensor components Ap given in Table IVl The values of Ap 
for 7th approximation fei^, 6or and for hard-disk molecules seems to be converged to three 
significant figures, as can be seen from Table Note that the off-diagonal components of 
Ap are zero, and Aff = — Ap is satisfied. Equation shows that the equation of state 
in the nonequihbrium steady state is not modified to first order. We indicate that the 
second-order pressure tensor should be uniform from the solubility conditions for the 
third-order solution 0^^^: 

•^-iVi ■ V/f ^dvi = 0. (56) 

Therefore, the pressure tensor Pij in eq. (|33j) becomes uniform since nnT is constant from 
the solubility conditions (jH)). 

We find that Ap for hard-disk molecules differs from that for the steady-state BGK 
equation and information theory not only quantitatively but also qualitatively: Aff < Ap^ 
for hard-disk molecules, while Aff = Ap' for the steady-state BGK equation and A|f > A^f for 
information theory. This kind of the difference has been also found in the three-dimensional 
case, jisl 

Each component of the kinetic temperature in the nonequihbrium steady state, i.e. 
for i = a; and y is calculated as 

which leads to 

T-^ = ni + >^P^l (58) 

for i = X and y. Values for the constants in the second-order term Ap are the same as Ap 
for i = j given in Table |3 For hard-disk molecules, Tj. becomes smaller than Ty regardless 
of the sign of Jx, which means that the motion of hard-disk molecules along the heat flux 
becomes dull. We note that T; for hard-disk molecules is isotropic to first order, that is, the 
equipartition law of energy holds. 
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The Shannon entropy in the nonequihbrium steady state 5* is defined via 



S = 



/ dc/log/ 



J — oo 



—riK log 



nm 



+ nK + Xs 




(59) 



l2TTKTi 



Values for the constant Xs are given in Table El Xs for 7th approximation bir, feor and b2r for 
hard-disk molecules seems to converge to four significant figures. It is found that Xs for hard- 
disk molecules is close to that for the steady-state BGK equation and information theory. 
This is because the second-order correction term in the Shannon entropy is determined only 
by the square of the first-order solution 0*^^^ where no important difference dependent on the 
kinetic equations or information theory appears. We note that the Shannon entropy in the 
nonequihbrium steady state is not modified to first order. 

VII. CONTRIBUTION OF THE STEADY HEAT CONDUCTION TO THE RATE 
OF CHEMICAL REACTION 

A. Calculation of the Rate of Chemical Reaction 

In the early stage of a chemical reaction between monatomic gas molecules: 



the rate of chemical reaction may not be affected by the existence of products, and the 
reverse reaction can be neglected. |62] From the viewpoint of kinetic collision theory, the 
chemical reaction rate (jUUI) can be described as 



for two dimension. Here f2 denotes the solid angle for two dimension. For the differential 
cross-section of chemical reaction (j{g), we have derived the line-of-centers model for the 
case of the two dimension. Its form becomes 



A + A ^ products. 



(60) 




(61) 




(62) 



with m mass of molecules and E* the threshold energy of a chemical reaction. 
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We calculate the rate of chemical reaction ()61|) with the two-dimensional line-of-centers 
model (jU^ using the explicit velocity distribution function of the two-dimensional steady- 
state Boltzmann equation for hard-core molecules to second order (|4Up . Substituting the 
expanded form of the velocity distribution function to second order in eq.(j21) into eq. (j61|) . 
we obtain 

i? = i?(°)+i?W + i?(2)^ (63) 
up to second order. The zeroth-order term of R, 

= JdvJdv,JdnJ f^'^ff^ga{g) = 2n'd ' e'^ , (64) 

corresponds to the rate of chemical reaction of the equilibrium theory. Similarly, the first- 
order term of R is obtained as 

= I rfv| rfv, I rfO I /W/f + <P?]ga{g), (65) 

where R^^^ does not appear because 0'-^^ is an odd functions of c, as is shown in eq. ()4Up . 
The second-order term of R, i.e. R^'^\ is divided into 

^^''^^ ^ I I ^""'1 '^^l (66) 



and 



R^'^'^ = Jdvldv,Jdnl + <t>?]gaig), (67) 

which exhibit the local nonequilibrium effect. Since the integrations (jMj) and (jU7j) have the 
cutoff form as in eq.()62|). the explicit forms of 0*^^-* and 0^^-' of the steady-state Boltzmann 
equation for hard-disk molecules are required to calculate i?*^^'^^ and R^'^'^\ respectively. 

In Fig. 121 we have confirmed that 0^^-* seems to converge to 7th Sonine approximation, 
so that we will show only the results calculated from 0*^^^ and 0*^^-' for 7th approximation 
of Sonine polynomials. In order to compare the results from the steady-state Boltzmann 
equation with those from the steady-state BGK equation and information theory, we also 
use the explicit forms of (p^^^ and 0*^^^ obtained in eq. ()46p for the steady-state BGK equation 
and in eq.()52|) for information theory. 

B. Local Nonequilibrium Effect on the Rate of Chemical Reaction 

Inserting 0*^^-' and 0*^^-' of eq.()40j) for the steady-state Boltzmann equation for hard-core 
molecules, eq.(|m|) for the steady-state BGK equation and eq.(|^ for information theory into 
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eqs. ()66p and ()67|) . and performing the integrations with the chemical reaction cross-section 
(jU^ . we finally obtain the local nonequilibrium effect on the rate of chemical reaction based 
on the line-of-centers model. The expressions of R^'^'^^ and R^"^'^^ become 

and 

respectively. The numerical values for and (3^ are listed in Tables IVTl and rVTTl respectively. 
As well as the three dimension case, the two-dimensional R^"^'^^ for steady-state Boltzmann 
equation is determined only by the terms of foor in ^'•^^ of eq. ()4U|) . This is because i?'-^'^-* 
of eq. (p7j) has x-y symmetry, so that the terms including (c^ — Cy) in (p^"^^ of eq. pUj) do not 
contribute to R^'^'^\ 

The graphical results of i?*^^-* compared with those of R^'^'^^ are provided in Figl^ Both 
of and i?(2,A) in FigElare scaled by n^/^dm^/^J^/K^/^T^/^. Note that jg the sum of 
R^"^'^^ and i?*^^'^) in eqs. ljfiHj) and ()fi9|) . As FigElshows, it is clear that R^"^'^^ plays an essential 
role for the evaluation of R^'^\ We have found that there are no qualitative differences among 
R^"^^ and R^'^'^^ for the steady-state Boltzmann equation, and those for the steady-state BGK 
equation and information theory, while they exhibit slight deviations from each other. The 
quantitative deviation in R^"^'^^ would not be observed if we adopted (p^^^ of the steady-state 
Boltzmann equation for the lowest Sonine approximation, because that is identical with the 
precise (p^^^ of the steady-state BGK equation and information theory. 




VIII. DISCUSSION AND CONCLUSION 



Fushiki has recently demonstrated that our analytical three-dimensional second-order 
solution of the steady-state Boltzmann equation for hard-core molecules agrees well with 
results of his numerical experiment u sing both a molecular dynamics simulation and a direct 
simulation monte carlo method, [ifll. l2c| Using the method developed in ref.[3], we have 
derived the velocity distribution function of the two-dimensional steady-state Boltzmann 
equation for hard-disk molecules explicitly to second order in the temperature gradient. 
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as was shown explicitly in eq. (pn|) and graphically in Figl^l We have calculated the two- 
dimensional equation of state for hard-disk molecules to second order from it. We believe that 
the second-order solution of the steady-state Boltzmann equation is physically important in 
that it reflects a nonequilibrium state far from equilibrium, called the local nonequilibrium 
state. 

We have found that there are qualitative differences between hard-disk molecules and the 
steady-state BGK equation in the nonlinear nonequilibrium transport phenomena based on 
the local nonequilibrium state: the second-order corrections appear for hard-disk molecules 
in the pressure tensor Pij and the kinetic temperature Ti, while no corrections to these 
quantities appear for the steady-state BGK equation, as Table |3 shows. This kind of 
qualitative differences was detected also in the three-dimensional case.^^ This discrepancy 
is due to the fact that (/-dependency cannot be absorbed in the single relaxation time of the 
BGK equation, which leads to the conclusion that the steady-state BGK equation neither 
capture the essence of hard-disk molecules nor possess the characteristics of any other models 
of molecules which interact with (/-dependency. We suggest that microscopic models which 
possess the property that its relaxation to the local equilibrium state is described only 
by a single relaxation time could not be applied to describe the nonlinear nonequilibrium 
transport phenomena. This suggestion may mean that the steady-state BGK equation could 
capture the essence of hard-disk molecules if one made the relaxation time depend on g or 
if one developed the steady-state BGK equation with multi-relaxation times. We note that 
the qualitative differences mentioned above still appear no matter which boundary condition 
is adopted, that is, the isotropy and the anisotropy of the pressure tensor in eq.(jSSI) and the 
kinetic temperature in eq. ()58p are not affected by any kinds of boundary conditions. 

We have examined information theory by the microscopic kinetic theory mentioned above, 
and consider the possibility of the existence of a nonequilibrium universal velocity distribu- 
tion function. The first-order velocity distribution functions for the steady-state Boltzmann 

is consistent with that 



equation for hard-disk molecules, i. e. the first-order terms in ecj. ^^^^^ 
derived by expanding Zubarev's velocity distribution ftmction jSClL l5lL l52|| . On the other hand, 
the explicit form of the second-order term in eq. (jlUI) definitely differs from the precise form 
for the steady-state BGK equation ()46|) or information theory ()52|). as FiglH shows. A 



information theory has been applied to nonequilibrium dilute gases j39|, |4c 



55, 



56, 



though 



57,^ 



we have found that information theory contradicts the microscopic kinetic models: all 
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the macroscopic quantities for information theory except for the Shannon entropy 5* in 
eq. ()59|) are quahtatively different from those for the steady-state Boltzmann equation and 
the steady-state BGK equation. These results indicate that characteristics of microscopic 
models appear in the local nonequilibrium state, that is, nonlinear nonequilibrium transport 
phenomena are sensitive to differences of kinetic models, so rather realistic models are needed 
when one investigates them. We can conclude that, though quite a few statistical physicists 



have believed the existence of a universal velocity distribution 



unction in the nonequilib- 



rium steady state by maximizing the Shannon-type entropv|lll. I47l l48l . l50l . l5lL l52l |53[, any 
universal nonlinear nonequilibrium velocity distribution function does not seem to exist in 
the two-dimensional well as the three-dimensional case, even when it is expressed 

only in terms of macroscopic quantities. We suggest that the entropy defined in eq. ljTTjl is 
not appropriate as the nonequilibrium entropy to second order though it is appropriate to 
first order, and that some nonequilibrium corrections dependent on microscopic models are 
needed for the nonequilibrium entropy to second order. 

The second-order solution of the steady-state Boltzmann equation for hard-disk molecules 
is indispensable for the calculation of the nonequilibrium effects on the rate of chemical 
reaction, since i?*-^^ does not appear and R^"^'^^ is remarkably larger than i?*^^'"^-* as FiglJ] 
shows. This indicates the significance of the second-order coefficients as terms which reflect 
the local nonequilibrium state. 
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APPENDIX A: CALCULATION OF Qkr 

From the definition of Qkr, ^kr can be calculated using the mathematical properties of 
the trigonometrical functions and Sonine polynomials. For example, cixQkr can be rewritten 

as 

f2KTY ^f^\^ Slicj) 2kT 



m 
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Integrating eq. ()Al|) over (2KT/m)i/2ci with /i from ea.dTD can be performed by using the 
following orthogonality properties. For Sonine polynomials, 

-l)P-«r(g + A; + 1) 



X'^e-^SU{X)Sl{X)dX 



for p = q and p = q + 1, and is zero otherwise. For the trigonometrical functions, 



(A2) 



2n 



COS n0 sm m 



0, 



(A3) 



and 



27r 



COS n0 COS m 



2n 



sm no sm m- 



(A4) 



'0 ^0 

with the Kronecker delta 6pq for m 7^ 0. We can calculate Dfc,r and E^^r defined as 

1- 



n n 

The results can be written as 



(A5) 



Dk,r = (4,0 + l)[{k + r + l)Bk+l^r - 5fc+l,r-l + (4,0 + l)5fc-l,r - (4,0 + l)(r + l)Bk-l,r+l] 



(A6) 



and 



-Efc,r = (4,0 + 1) [{k + r + l)Ck+l,r — Ck+l,r-l — Ck-l,r + (r + l)Cfc_i,r+l] 



(AT) 



Additionally, Ci^d^Qkr can be rewritten as 



kr 



T 



<5fc,r-l(Cl) - (r + -)Qk,r 



(A8) 



Therefore, by integrating eq. ()A8|) over (2KT/m)^/^Ci, with fi from eq. ()14|) . it is found that 



C^x^xQ kr ^ 



dxT 
T 



Dk,r-i - {r + -)Dk,r 



(A9) 



Similarly CiydyQkr is obtained by replacing the differential coefficients with respect to x by 
the corresponding differential coefficients with respect to ?/, the D^^s by the corresponding 
Ek^r, respectively. Substituting these results into eq.(j2Sl), ^kr finally becomes eq.(|^. 
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APPENDIX B: CALCULATION OF F^V(x) 
1. Details of 

The details of F^^{x) are written in this Appendix. Substituting the general forms of /i, 
/2 in eq. (fT^ and Q'^^ in eq. ipTj) into Fl^{x) in eq. lj^ . -Ffcr(x) can be written as 

ni,n2,fci,fc2 

where S^'^;^"^^"^ (x) is the characteristic integral defined as 

sS'^rW = / |exp[-(c? + c^)]n(c;)n,(cOn,(c2)5I(cf)5-(c?)5-(c^)^dc2dc^ 

(B2) 

and 

(B3) 

We note that values of the latter is obtained from those of the former by a transformation 
of axes, and that '^k'okl'^^ix) = '^k'k^o"^^ix) = from eq. ffTHj) . The integral containing 
Yfei(ci)2'fc2(c2) + Zfc^(ci)Yfc2(c2) becomes zero, owing to the symmetry of the trigonometrical 
functions. The factor W^'^^^'J!^ in eq. ()Bl|) is defined as 

An"^ / m ^ — ' — - 



which is obtained from the prefactors and the coefficients in the general form of /i, /2 in 
eq. (P|) and Q^^ in eq.^. 

We find that it is necessary only to evaluate the characteristic integral ■^/-'^.^ 
order to calculate F^j.{x)- Our calculation of S^'^'"^^"^ (x) is written in the next subsection. 
Once the characteristic integral '^'k'k^k^^ ix) has been derived, F^^{x) is now calculated from 
eq.pq) with l^fcXX' in eq.(|B4)). 

2. Calculation of 

We shall explain how to calculate '^^'k^k-^^ix) which appears in ea. ()B2|) . The calculation 
has been performed mainly based on the method developed in ref. 
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a. Introduction of Qk MM 



Using eq.(fl3j). the characteristic integral corresponds to the coefficient of 

©fc,fci,fe2 = i^kMM j J Yk{c[)Yk^{ci)Yk2{c2)exp{-{acl + Pcl + -fcf)}gdc2dci, (B5) 



that is 



0,,,,,,^ J2 ^l:^X\x)s^t-^u-^. (B6) 



r,ni,n2 



In eq. ()B5|) a, /? and 7 are defined as 



« = = 7=T^, (B7) 

1 — t \ — u 1 — s 

and z/fc,fci,fc2 is given by 



i^fcA,;^^ = (1 - - ty-\l - u)-'^-' , (B8) 

\ m J 

and 

Yfc(ci) = cf cosA;0, rfci(ci) = c^^ cos and ^^^(cs) = 4' cos A;202- (B9) 
Finally, we need only to evaluate the characteristic integral QkMM order to calculate 

■~Y,r,ni,n2 / \ 
^k,ki,k2 \X,)- 

b. Derivation of the Inductive Equation 



In order to evaluate QkMM eq. ()B5|) . let us derive an inductive equation for QkMMi 
which is related to QkMM by 

^k,kiM = ^kMM^kMM- (BIO) 

By replacing Ci^ and C2x by Ci^ — w and C2x — ui, respectively, and will be changed to 
^ix ^2x At the same time, the relative speed g is not modified, and the value of 

Qk,kiM is unchanged. Therefore, QkMM is independent of w and differentiation of QkMM 
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with respect to w gives zero. After this differentiation has been performed and w has been 
set to zero, it is found that 




exp{ — (ac^ + [3c^ + '^c''^)}gdc2dci x 
janYfe^ (^fei+i + clYk,-i){5k,fl + 1) + pYkYk,iYk.,+i + c^n^-Ol^a.o + 1) 

(Bll) 

by using the formulae 

Ci.n,(ci) = {n,+i(cO + cfn,_i(cO} , (B12) 

and 

^^^ = hYk.Mci), (B13) 
for i = 1 and 2. From eqs. ()B5|) and ()B10|) . eq. ()Bll|) leads to the inductive equation 

a(4i,0 + l)0fc,A:i+l,A:2 " "(4i,0 + 1) " ^'^^ 

fc A; 1 

+ /3(5fe2,0 + l)0fc,fci,fc2 + l - /5(5fc2,0 + 1) 

+ 7(4,0 + l)efc+i,.,.2 - 7(4.0 + 1)^^^^^ 

c'7 



(B14) 



for Qk,kiM- Because of this inductive equation, once the initial value Qo^kiM is known for 
all ki and k2, then the values of the integral Qk,kiM for k, ki and k2 can be obtained, 
and QkMM is then obtained from eq. ()B10|) . 

c. Calculation of the Initial Value 

In principle, the initial value of the inductive equation ()B14|) . Qo^kiMi ^^"^ ^le obtained 
and written explicitly from eq. ()B5|l . changing the variables Ci and C2 to V = (ci + C2)/2 
and g = Ci — C2. Though we have directly calculated Oo,fci,A:2 oiily fo^' ki = k2 = 0, 1, 2, 3, 
4, 5, 6, 7 and 8, they are sufficient to get all the results shown in Appendices O and |D| We 
do not show the explicit expressions of the initial values in this paper because they are too 
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complicated. We have also confirmed that the initial value Qqmm becomes zero for ki ^ k2. 
Note that Qqmm is obtained by i^omm^omm ^o"^ ea.()Bl();i. 



d. Evaluation of E,^'^'^l'^'^{x) 

Using the inductive equation ()B14|1 and the initial value Qqmm calculated in Sec JB 2 c( 
the values of the integral QkMM ^2, can be obtained with the relation 

()Bin|l . The result is that QkMM vanishes k = |fci — A;2| + 2g, where g is a positive integer or 
zero. In order to obtain it is sufficient to have QkMM o^^Y ^1 — ^2- This 

is because the value of '^'k^i'^'^ix)^ with ki and ^2, rii and n2 interchanged, corresponds to 
the value of S^'^:^"^^"^ (x) with x replaced by vr - x- Thus, if ki ^ k2 or m ^ n2, '^Ifk^'k^^x) 
for any set of ki and /c2, ni and n2 corresponds to 

-kMM yx)=^kMM yx)+^kMM " x), (B15) 
for ki > k2 ; if ki = k2 and rii = 77,2, then '^^'I'^k-^'^ (x) gives the required value at once. 



APPENDIX C: CALCULATION OF THE FIRST ORDER COEFFICIENTS 

Let us explain how to obtain the first-order coefficients, that is, how to solve the integral 
equation (jS)). To begin with, we calculate fl^^ in eq.(j2Sl) to first order; for first order 
corresponds to the right-hand side of eq.©. It can be calculated only by substituting 
Bqq = 1 into the expressions of Dk^r and Ek^r in eqs. ()A6|) and ()A7|) : the coefficient Bqq = 1 
corresponds to /i = fi'^\ and no higher-order terms appear in to first order. It finally 
becomes 

Now fl^j. for first order is found to vanish unless = 1, so that we need calculate only 
A^^ for first order; as was mentioned in the end of Seel ^1 we do not need to consider the 
case in which the right-hand side of eq.(jS} becomes zero. 2] To derive A^^ in eq. ()28|l for first 
order, we must calculate both W^'^^'J!'^ and S^'^J'^^ in F^^^^x) of eq.(jBT} for first order, as 
was shown in Appendices IB II and IB 21 The result for A^^ to first order can be written finally 
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in the form 

where the set of the coefficients B\^^Bqq is obtained from W['^^ff in eq. ()B4|) . To first order, 
/i in eq.(|^ contains only Bqq = 1 and the first-order coefficients Bl^^_^ and C^^„^; /2 in 
eq. lf^ also contains Bqq = 1, Bl^^^ and Cl^^^ to first order. Thus, we obtain only the term 
Bl^^BoQ from WI'^q^ to first order using the fact that Flj.{x) = unless k = \ki — + 2q. 
Note that it is sufficient to consider only the case for ki > k2 as is explained in Appendix 
IB 2\ and that we set g = 0. The matrix M^f'J^^'^ is thus obtained 

^1,1,0 = n-^Tr{r + 2) l [^^'1'° W-"M-o (0)]sm-dx, (C3) 
using eqs.(I2Hl), jHH) and (fR4jl . 

For /c = 1, eq. ()H()|l gives a simultaneous equation determining the first-order coefficients 

ni>l 

from eqs. ()Cl|) and ()C2|) . Note that we need only to obtain the first-order coefficients for 
ni > 1, because -Bio = from eq.lfT^. We have calculated the matrix M^^'J'^'^ for 1 < r < 7 
and 1 < TT-i < 7 from eq. ()(y3jl . and we have also confirmed that M^f'Q^'^ for 1 < rii < 7 
calculated from eq.((nHl) vanishes. Our result for M^f'J^^'^ for 1 < r < 7 and 1 < < 7 is 
given in Appendix |E1 At last, we can determine the first-order coefficients B\^_^ by solving 
the simultaneous equation ()C4|) . that is, can be obtained as 

BU = J2^l{Ml',y'')-\ (C5) 

where represents the inverse matrix of a matrix X. Finally, the results of the first-order 
coefficients -B^^^^, i-e. the first-order in eq. lfT^ . can be calculated as in eq.dHU- 



APPENDIX D: CALCULATION OF THE SECOND ORDER COEFFICIENTS 

rII 

^kr 

We explain how to obtain the second-order coefficients, that is, how to solve the integral 
equation ®. The coefficients of first order, i.e. B],^ and C\^, have been obtained as are 
given in eq.(|HT|. so that we can employ them to determine the second-order coefficients. 

27 



To begin with, we calculate fif^ in eq.(j2SI) for second order; fl^^ for second order corre- 
sponds to the first term on the right-hand side of eq.®- It can be calculated by substituting 
Bl^. and Cl^ into the expressions of D^^r and E^^r in eqs. fjA6|) and ()A7|) : no other terms ap- 
pear in Q^j. for second order. The results of the tedious calculation of to second order 
finally become as follows. For k = 0, Q^j. becomes 



nl = 0, 



(Dl) 



for r = and 1, 



27r2 /2kT\2 



m 



{VT)'{^h2-lbn)+T{V'T) {3b 



12 



'II, 



for r = 2, and 



2ti^ /2kT 



{T(V'T) [(r + - 



+ {VTf 



:d2) 



(D3) 



for r > 3. Note that values for the constants hir are summarized in Table HI For k = 2, Qfj. 
becomes 

/2/tTV r 6ii 



20 rfr2 V m 



(D4) 



and 



and 



-fell - 36i2 



+ T[dlT-d^T] [611-2612] 



(D5) 



2r 



TT 



2kT\ 2 



2 



+ T [dlT - 9jT] [-(r + l)6i,,+i + 61,]} 



1 1 
-(r + -)(r + l)6i,,+i + (2r + -)6i, - 6i,,_i 



(D6) 



for r > 2. For = 1 and > 3, we find that fi^^ for second order becomes 



C = 0, 



(D7) 
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for any value of r. 

Next let us calculate Af^ in eq.(j2HI) for second order. In order to derive Af^ for second 
order, we have to calculate Wl']^^'J!^ and '^^'l'^^^^ in F^^{x) of eq. ()Bl|) to second order, as 
was shown in Appendix IB 11 For k = 0, A^^ to second order results in 

K = Boo E Bi\Mi-;--'+ J2 <,<.Mo^i:r'- + cucUKiT'""'- m 

ni>2 "1,^2 ^1,1^2 

B\,i^_^ from fi and B\^^ from /2 of the set of the coefficients Bl^^B\^^ are the first-order 
coefficients obtained in eq. (jHT|) . so that Bl^^B\^^ is second order. Similarly, C\^^C\^^ is also 
second order. The second and the third terms on the right-hand side of eq. flDSj) correspond 
to J(/i, /2) in the integral equation (jHl). To second order, /i of eq. lj^ contains only Boo = 
1, and C^.„. obtained in eq.(jSH), and -B"„. to be determined here for i = 1 and 2. 

Therefore, we can only obtain the sets of the terms in eq. ()D8|) for second order by using 
the fact that Fl^{x) = unless k = — + 2g. We should derive the second-order 
coefficients B^^_^ only for rii > 2, because Boo = 1 and Boi = from eq.lfT^. Note that it 
is sufficient to consider only the case for ki > k2, as is explained in Appendix IB 2\ and that 
-^LjC'inj + CiniBin2 does uot appear, (see Appendix IB Ij) 
The matrix Mq '^^^q^'^ in eq. ()D8|l is obtained as 

MY,r,ni,0 ^u/t / i^Y,r,ni,0/ \ ^Y,r,ni,0 /r,\l ■ A i /-nvnN 

0,0,0 = / Fo,b;o (x) - ^0,0,0 (0)] -dx, (D9) 

using eqs. (PH|) . ()B1|) and ()B4|) . Similarly, the matrices M^^'"^'"^ in eq. ()D8|l are derived as 
Mo^Cf" = [2£'r'"^(x)-Sj[;r'-(0)]sin|dx, 

(DIO) 

while we have confirmed Mo^l'f^'"' = Mo^i'f^'"^ Equations (|D8l), ^ and (|D3l) lead 

to a simultaneous equation to determine the second-order coefficients -BgJij : 

= E 1^0. - E (^in.^^n. + C^L.C^LJ (Mj---'^)-. 

(Dll) 

We have calculated the matrix Mqq'q^''^ for 2 < r < 6 and 2 < < 6 from eq. ()D9|) . and also 
confirmed that Mqq'q^'^ vanishes for r = 0, 1 and 2 < ni < 6 or for 2 < r < 6 and rii = 0, 1. 
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We have calculated the matrix Md^f'"^'"^ for2<r<6, 1 < rii < 7 and 1 < n2 < 7 from 
eq. ()D10|) . and confirmed M|}'f'"^'"^ vanishes for r = 0, 1, 1 < ni < 7 and 1 < n2 < 7. Our 
results for mI'q^^^'° for 2 < r < 6 and 2 < rii < 6 and Mg^J^f for 2 < r < 6, 1 < ni < 7 
and 1 < "^-2 < 7 are given in Appendix [Ej Finally, we can determine the second-order 
coefficients B^^_^ in /i, i.e. the second-order Bqj. in eq. lfT^ as in eq.(jSni)- 
Similarly, for k = 2, A^^ for second order results in 

ni>0 ni,«2 iii,n2 

using the fact that F^^{x) = unless k = \ki — + 2q. Note that we have confirmed 
M^q'q^'^ becomes zero. The matrix mJ^^'^'J'^'" in eq. ()D12jl is obtained as 

^2,2,0 = fe.o (x) - ^2:2,0 (0)] sm -dx, (D13) 

using eqs. ipHj) . ()B1|) and ()B4|) . Similarly, the matrices M^^'^^^'"^ in eq. ()D12|) are derived as 

^™2 2 /•27r 

^cr'"' = / [sar^'"^(x)-sar^'"^(o)]sin|dx, 

•/ 

(D14) 

while we have confirmed M^'f""' = -M^;l'^^''^\ Thus, eqs. flDT2|l . (iDill . (lD5|l and (jD6)l 
lead to a simultaneous equation to determine the second-order coefficients B^^^: 

In order to derive the second-order coefficients B^^^ for rii > 0, we have calculated the matrix 
M2 2'o^'^ for < r < 6 and < rii < 6 from eq. ()D13|) . and also the matrices M^f'"^'"^ for 
0<r<6,l<rii<7 and 1 < ri2 < 7 from eq. ()D14|l . Those results are given in Appendix 
IeI The second-order coefficients Bj}^^_^ in /i, i.e. the second-order B^j. in eq.(fT^. can be 
written in the final form shown in eq.()37|). 

We need to consider eq. ()30|) only for k = and 2 for second order: it is not necessary 
to consider ea. (l3UI) f or even k furthermore, which was first expected in ref. and recently 
confirmed in ref.[l9|. For odd fc, n^^, to second order is found to be zero, and no terms 
corresponding to J(/i, /2) in the integral equation ©, i.e. the second and the third terms 
on the right-hand side of eqs. ()D11|) or ()D15|) appear, so that any second-order terms Bj}^ do 
not appear for odd /c.jist 
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APPENDIX E: MATRIX ELEMENTS 



1. M^'^^Q^'^ for 1 < r < 7 and 1 < ni < 7 

The matrix elements M^'['q^'^ for 1 < r < 7 and 1 < "n-i < 7 divided by M^^^q'* 
-2dn^ (^^) ^ calculated from eq. ()C3|l are given as follows. 

/ 1 -8.333 X 10-2 -2.604 x 10"^ -1.302 x lO""* -6.782 x 10-" -3.391 x 10"^ -1.589 x lO"* \ 

-5.000 X 10^1 1.625 -1.185 x 10"^ -3.451 x 10"^ -1.648 x 10""' -8.308 x lO"'^ -4.053 x 10"^ 

-1.875 X 10^1 -1.422 2.165 -1.327 x 10"^ -3.466 x 10"^ -1.531 x lO""* -7.279 x 10"^ 

-1.875 X 10-1 -8.281 x 10"! -2.654 2.645 -1.383 x lO'^ -3.234 x 10"^ -1.314 x lO-* 

-2.930 X 10^1 -1.187 -2.079 -4.148 3.081 -1.398 x 10"^ -2.94716687 x 10"^ 

-6.152 X 10^1 -2.512 -3.859 -4.074 -5.873 3.483 -1.395 x lO'^ 

V -1.615 -6.864 -1.027 x 10 -9.273 -6.932 -7.810 3.858 / 



2. M, 



Y,r,ni,0 



for 2 < r < 6 and 2 < m < 6, and M^J^f ^'^^ for 2 < r < 6, 1 < ni < 7 and 



0,0,0 



1 < n2 < 7 



The matrix elements M^'^^'q^'^ for 2 < r < 6 and 2 < ni < 6 divided by M^qq'^ 
-Adn^ (^^) ^ calculated from eq. ()D9|) are given as 



1 

-7.500 X 10-1 
-3.750 X 10-1 
-4.688 X 10-1 
-8.789 X 10-1 



-8.333 x 10-2 -2.604 x 10-3 -1.302 x 10-"^ -6.782 x 10 

1.688 -1.152 x 10-1 -3.223 x 10-3 -1.495 x 10-^ 

-1.844 2.251 -1.287 x 10-1 -3.204 x 10-3 

-1.289 -3.217 2.741 -1.342 x 10-1 

-2.153 -2.883 -4.833 3.182 

fY,2,ni,n2 



y,2,l,l 



The matrix elements Mq'^'^^'""^ for 1 < ni < 7 and 1 < n2 < 7 divided by Mq[\ 
■ ^Jrj^i (^^) ^ calculated from eq. ()D10|) are given as 



/ 1 5.000 X 10-1 0.000 

5.000 X 10-1 3.750 x 10"! 9.375 x 10"! 

0.000 9.375 X 10-1 1.230 

-9.375 X 10-1 8.203 x lO-i 5.537 

-4.102 -1.846 1.015 x 10 

-1.661 X 10 -2.030 X lOi 0.000 

V -7.106 X 10 -1.320 X 10^ -1.856 x 10^ 



-4.102 -1.661 X 10 -7.106 x 10 \ 

-1.846 -2.030 X 10 -1.320 x 10^ 

1.015 X 10 0.000 -1.856 x 10^ 

6.598 X 10 1.856 x 10^ 2.629 x 10^ 

1.547 X 10^ 1.315 X 10^ 4.995 x 10^ 

1.315 X 10^ 3.746 x 10^ 3.934 x lO^* 

2.629 X 10^ 4.995 x 10^ 3.934 x lO"* 1.319 x lO'^ J 



-9.375 X 10-1 
8.203 X 10-1 
5.537 
1.015 X 10 
6.598 X 10 
1.856 X 10^ 



31 



The matrix elements M(^J"i'"' for 1 < m < 7 and 1 < ria < 7 divided by M^^'^^l'^ are 
given as 





-2.500 X 10-1 


1.125 


5.000 X 10-1 


-4.688 X 10-2 


-2.490 


-1.482 X 10 


-7.913 X 10\ 




1.125 


2.813 X 10-1 


7.969 X 10-1 


1.436 


2.000 


-4.153 


-7.360 X 10 




5.000 X 10-1 


7.969 X 10-1 


7.471 X 10-1 


3.845 


1.130 X lOi 


3.045 X lOi 


2.990 X lOi 




-4.688 X 10-2 


1.436 


3.845 


5.768 


4.187 X 10 


1.681 X 10^ 


6.534 X 10^ 




-2.490 


2.000 


1.130 X 10 


4.187 X 10 


8.506 X 10 


7.926 X 10^ 


4.009 X 10^ 




-1.482 X 10 


-4.153 


3.045 X 10 


1.681 X 10^ 


7.926 X 10^ 


2.021 X 10^ 


2.295 X 10-* 


v 


-7.913 X 10 


-7.360 X 10 


2.990 X 10 


6.534 X 10^ 


4.009 X 103 


2.295 X 10* 


7.027 X 10* y 



The matrix elements Mo^^'^f ^'^^ for 1 < m < 7 and 1 < ria < 7 divided by M^^l^l'^ are 
given as 



/ -1.302 X 10-2 
-2.305 X 10-1 
7.813 X 10-1 
4.946 X 10-1 
-8.331 X 10-2 
-4.567 
V -3.443 X 10 



-2.305 X 10-1 
2.295 X 10-1 
3.960 X 10-1 
7.784 X 10-1 
1.869 
3.623 
-6.020 



7.813 X 10-1 4.946 x lO-i 
3.960 X 10-1 7.784 x lO-i 



3.062 X 10-1 
1.513 
4.918 
1.817 X 10 
6.396 X lOi 



1.513 
2.111 

1.527 X 10 
6.642 X 10 
3.216 X 102 



-8.331 X 10-- 
1.869 
4.918 
1.527 X 10 
2.971 X in 
2.782 X 102 
1.502 X 10^ 



-4.567 
3.623 
1.817 X 10 



-3.443 X 10\ 

-6.020 
6.396 X 10 



6.642 X 10 3.216 x 102 

2.782 x 102 1.502 x 10=^ 

6.882 X 102 7.871 x 10^ 

7.871 X 10^ 2.355 x 10* / 



The matrix elements Mo^f]"i'"' for 1 < rii < 7 and 1 < ria < 7 divided by 'j^^"'" are 
given as 



r^,2,i,i 



/ -9.115 X 10-4 -1.074 X 10-3 -1.289 x lO-i 5.948 x lO-i 4.914 x lO-i 



-1.074 x 10^2 
-1.289 x 10-1 
5.948 X 10-1 
4.914 X 10-1 
-1.037 X 10-1 
\ -7.105 



-4.272 X 10^2 
2.509 X 10-1 
3.006 X 10-1 
7.575 X 10-1 

2.279 

5.713 



2.509 X IQ-i 
1.104 X 10-1 
4.718 X 10-1 

1.484 

5.886 
2.621 X 10 



3.006 X 10-1 7.575 x lO-i 
4.718 x 10-1 1.484 



6.042 X 10-1 
4.226 
1.845 X 10 
9.461 X 102 



4.226 

7.918 
7.326 X 10 
4.000 X 10 



-1.037 X 10-1 
2.279 
5.886 
1.845 
7.326 X 10 
1.772 X 102 



-7.105 \ 

5.713 
2.621 X 10 
9.461 X 10 
4.000 X 102 
2.017 X 10 



The matrix elements M(^J|^f^^'"^ 
given as 



2.017 X 10^ 5.941 X 10^ / 

y,2,i,i 



for 1 < ni < 7 and 1 < ^2 < 7 divided by Mq ^ ' are 





-6.104 X 10 


-5 


-6.978 X 10-4 


-5.046 X 10-3 


-8.158 X 10-2 


4.782 X 10-1 


4.898 X 10-1 


-1.087 X 10- 


i\ 




-6.978 X 10 


-4 


-1.976 X 10-3 


-3.976 X 10-2 


1.575 X 10-1 


2.407 X 10-1 


7.397 X 10-1 


2.678 






-5.046 X 10 


-3 


-3.976 X 10-2 


5.664 X 10-2 


1.387 X 10-1 


3.788 X 10-1 


1.450 


6.802 






-8.158 X 10 


-2 


1.575 X 10-1 


1.387 X 10-1 


1.505 X 10-1 


9.806 X 10-1 


4.159 


2.137 X 10 






4.782 X 10- 


1 


2.407 X IQ-i 


3.788 X 10-1 


9.806 X 10-1 


1.749 


1.579 X 10 


8.557 X 10 






4.898 X 10- 


1 


7.397 X 10-1 


1.450 


4.159 


1.579 X 10 


3.727 X 10 


4.194 X 102 




v 


-1.087 X 10 


-1 


2.678 


6.802 


2.137 X 10 


8.557 X 10 


4.194 X 102 


1.217 X 103 


/ 
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3. M2;2fl''^ for < r < 6 and < ni < 6, and M^^l'^^'""^ for < r < 6, 1 < ni < 7 and 

1 < na < 7 

The matrix elements mJ^'^'J^^'^ for < r < 6 and < ni < 6 divided by M2 2q'^ = 
—Adin? (^^) ^ calculated from eq. ()D13|) are given as 



/ 




1 




-S.333 X 10- 


2 


-2.604 X 10-» 


-1.302 X 10" 


4 


-6.7S2 X 10" 


6 


-3.391 X 10- 


7 


-1.5S9 X 10- 


8\ 




-2 


500 X 10- 


1 


1.063 




-S.OOS X 10-2 


-2.376 X 10" 


3 


-1.149 X 10" 


4 


-5.S49 X 10" 


6 


-2.S74 X 10" 


7 




-6 


250 X lo- 


2 


-6.406 X 10- 


1 


1.232 


-7.939 X 10- 


2 


-2.143 X 10" 


3 


-9.707 X 10- 


5 


-4.706 X 10- 


6 




-4 


ess X 10" 


2 


~2.S52 X 10" 


1 


-1.191 


1.420 




-7.S25 X 10" 


2 


-1.900 X 10- 


3 


-7.953 X 10" 


5 




-5 


S59 X 10" 


2 


-3.311 X 10- 


1 


-7.713 X 10-1 


-1.S7S 




1.607 




-7.676 X 10- 


2 


-1.6S0 X 10" 


3 




-1 


025 X 10" 


1 


-5.S96 X 10" 


1 


-1.223 


-1.596 




-2.6S7 




1.7S7 




-7.510 X 10" 


2 


\ 


-2 


307 X 10" 


1 


-1.391 




-2.S46 


-3.207 




-2.S22 




-3.605 




1.961 





The matrix elements M^f^"^'"^ for 1 < "n-i < 7 and 1 < ^2 < 7 divided by M^fi'^ 
i^)^ {{dxTf - {dyTY} calculated from eq. (fm4|l are given as 



16dT2 



/ 1 


1.500 


1.875 


3.2S1 


7.3S3 


2.030 X 10 


6.59S X 10 \ 


1.500 


9.375 X 10-1 


3.2S1 


7.3S3 


2.030 X 10 


6.59S X 10 


2.474 X 10^ 


1.S75 


3.2S1 


3.691 


2.030 X 10 


6.59S X 10 


2.474 X 10^ 


1.052 X 10^ 


3.2S1 


7.3S3 


2.030 X IQi 


3.299 X IQi 


2.474 X 10^ 


1.052 X 10^ 


4.995 X 10=* 


7.3S3 


2.030 X 10 


6.59S X 10 


2.474 X 10^ 


5.25S X 10^ 


4.995 X 10^ 


2.622 X 10"' 


2.030 X 10 


6.59S X 10 


2.474 X 10^ 


1.052 X 10^ 


4.995 X 10^ 


1.311 X 10^ 


1.50S X 10^ 


V6.59S X 10 


2.474 X 10^ 


1.052 X 10^ 


4.995 X 10^ 


2.622 X 10"' 


1.50S X 10^ 


4.712 X 10^ J 



The matrix elements M2 { [ 
nven as 



Y,l,ni,n2 



y,o,i,i 



for 1 < ni < 7 and 1 < ^2 < 7 divided by [ \ 



are 



/ 2.250 


2.125 


3.094 


6.680 


1.S25 X 10 


5.952 X 10 


2.246 X 


102\ 


2.125 


7.969 X 10-1 


2.930 


8.40S 


2.999 X 10 


1.231 X 10^ 


5.650 X 


10^ 


3.094 


2.930 


2.563 


1.523 X 10 


6.21s X 10 


3.011 X 10^ 


1.624 X 


10^ 


6.6S0 


S.40S 


15.23 X 10 


2.094 X 10 


1.691 X 10^ 


S.815 X 10^ 


5.324 X 


10^ 


1.S25 X 10 


2.999 X 10 


6.21s X 10 


1.691 X 10^ 


3.170 X 102 


3.221 X 10=* 


2.029 X 


10=* 


5.952 X IQi 


1.231 X 10^ 


3.011 X 10^ 


S.815 X 10^ 


3.221 X 10^ 


7.649 X 10^ 


9.343 X 


10* 


\ 2.246 X 10^ 


5.650 X 10^ 


1.624 X 10^ 


5.324 X 10^ 


2.029 X lO"! 


9.343 X 10'' 


2.6S6 X 


lO^y 



The matrix elements M2 '{ i 
nven as 



Y,2,ni,n2 



y,o,i,i 



for 1 < Til < 7 and 1 < n2 < 7 divided by M2 \ [ 



are 
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/ -3.211 


3.855 


2.370 


4.590 


1.368 X 10 


5.090 X 10 


2.199 X 10^^ 


3.855 


8.413 X 10-1 


2.059 


4.478 


1.522 X 10 


6.873 X 10 


3.660 X 10^ 


2.370 


2.059 


1.799 


8.862 


2.974 X 10 


1.351 X 10^ 


7.762 X 10^ 


4.590 


4.478 


8.862 


1.268 X 10 


9.068 X 10 


4.030 X 10^ 


2.253 X 10^ 


1.368 X 10 


1.522 X 10 


2.974 X 10 


9.068 X 10 


1.778 X 10^ 


1.650 X 10^ 


9.068 X 10^ 


5.090 X 10 


6.873 X 10 


1.351 X 10^ 


4.030 X 10^ 


1.650 X 10^ 


4.090 X 10^ 


4.651 X 10"' 


2.199 X 10^ 


3.660 X 10^ 


7.762 X 10^ 


2.253 X 103 


9.068 X 103 


4.651 X 10" 


1.391 X lO^ ) 



The matrix elements M^f '^''""^ for 1 < rii < 7 and 1 < ria < 7 divided by M^f 'l'^ are 
given as 



/ 2.493 X 10^1 


-3.585 


3.119 


2.456 


5.985 


2.271 X 10 


1.070 X 10^ 


\ 


-3.585 


1.265 


1.474 


2.290 


5.928 


2.340 X 10 


1.257 X 10^ 




3.119 


1.474 


9.898 X 10-1 


4.274 


1.257 X 10 


4.831 X 10 


2.443 X 10^ 




2.456 


2.290 


4.274 


5.626 


3.800 X 10 


1.565 X 10^ 


7.774 X 10^ 




5.985 


5.928 


1.257 X 10 


3.800 X 10 


7.156 X 10 


6.429 X 10^ 


3.357 X 10=* 




2.271 X 10 


2.340 X 10 


4.831 X 10 


1.565 X 10^ 


6.429 X 10^ 


1.557 X 10^ 


1.733 X 10" 




1.070 X 10^ 


1.257 X 10^ 


2.443 X 10^ 


7.774 X 10^ 


3.357 X 10^ 


1.733 X 10" 


5.108 X 10" 


/ 



The matrix elements M211 
given as 



Y,4,ni,n2 
,1,1 



for 1 < ni < 7 and 1 < 712 < 7 divided by M^fi'^ are 



/ 7.719 X 10"^ 2.301 X 10" 



2.301 X 10-1 
-2.277 
2.579 
2.484 
7.307 
\ 3.369 X 10 



-8.267 X 10-1 
1.678 
1.262 
2.404 

7.2919 
3.285 X 10 



2.579 2.484 7.307 3.369 x 10 \ 

1.262 2.404 7.292 3.285 x 10 

1.778 4.645 1.607 x 10 7.049 x 10 

2.037 1.276 X 10 5.015 x 10 2.363 x 10^ 

1.276 X 10 2.262 x 10 1.960 x 10^ 9.990 x 10^ 

5.015 X 10 1.960 X 10^ 4.589 x 10^ 5.000 x 10^ 



-2.277 
1.678 
5.255 x 10-1 
1.778 
4.645 
1.607 X 10 

7.049 X 10 2.363 x 10^ 9.990 x 10^ 5.000 x 10^ 1.442 x lO" / 



The matrix elements M^^'^'"''^ for 1 < ni < 7 and 1 < < 7 divided by Mlf^^'^ are 
given as 



/ 3.858 X 10" 


4 


6.369 X 10-=* 


1.235 X 10^1 


-1.567 


2.185 


2.488 


8.576 \ 


6.369 X 10" 


3 


4.532 X 10^2 


-8.947 X 10^1 


1.178 


1.088 


2.463 


8.595 


1.235 X 10" 


1 


-8.947 X 10^1 


4.736 X 10-1 


7.623 X 10-1 


1.573 


4.855 


1.940 X 10 


-1.567 




1.178 


7.623 X 10-1 


6.743 X 10-1 


3.723 


1.366 X 10 


6.171 X 10 


2.185 




1.088 


1.573 


3.723 


6.063 


4.991 X 10 


2.467 X 102 


2.488 




2.463 


4.855 


1.366 X 10 


4.991 X 10 


1.119 X 10^ 


1.186 X 10^ 


V 8.576 




8.595 


1.940 X 10 


6.171 X 10 


2.467 X 10^ 


1.186 X 10^ 


3.331 X 10^ / 



The matrix elements M^^f f"^'"^ for 1 < ni < 7 and 1 < n2 < 7 divided by M^f^^'^ are 
given as 
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/ 2.012 X 10-5 2.953 x 10"* 



2.953 X 10- 


4 


1.169 X 10- 


3 


2.997 X IQ- 


3 


4.239 X IQ- 


2 


7.361 X 10" 


2 


-5.364 X 10 


-1 


-1.141 




8.671 X 10- 


1 


1.888 




9.507 X 10- 


1 


2.482 




2.494 





2.997 X IQ- 


3 


7.361 X 10-2 


-1.141 


1.888 


2.482 


4.239 X 10- 


2 


-5.364 X 10-1 


8.671 X 10-1 


9.507 X 10-1 


2.494 


-2.109 X 10 


-1 


5.791 X 10-1 


5.711 X 10-1 


1.395 


4.977 


5.791 X 10- 


1 


2.404 X 10-1 


1.035 


3.364 


1.422 X 10 


5.711 X 10- 


1 


1.035 


1.473 


1.122 X 10 


5.291 X 10 


1.395 




3.364 


1.122 X 10 


2.374 X 10 


2.425 X 10^ 


4.977 




1.422 X 10 


5.291 X 10 


2.425 X 10^ 


6.591 X 10^ 
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FIG. 1: Schematic description of an interaction in two dimension. 
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FIG. 2: The scaled ^*^^)s for hard-disk molecules. The dash-dotted line, the dotted line, the short- 
dashed line, the long-dashed line and the solid line correspond to the scaled cj)^'^^ for hard-disk 
molecules with the 3th, 4th, 5th, 6th and 7th approximation 6or and b2r, respectively. Note that 
we put Cy = 0. 
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FIG. 4: Direct comparison of the scaled <^^^^ for hard-disk molecules with those for the steady- 
state BGK equation and information theory. The solid line, the dashed line and the dash-dotted 
line correspond to the scaled ^^^^ for hard-disk molecules with 7th approximation 6or and h2ri the 
steady-state BGK equation and information theory, respectively. Note that we put = 0. 
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FIG. 5: Scaled R^"^^ compared to scaled R^"^'^^ as a function of E* / kT for the two-dimensional 
case. The solid line, the long-dashed line and the dotted line show R^"^^ for hard-disk molecules, 
the steady-state BGK equation and information theory, respectively. The dashed line and the dash- 
dotted line represent i?^^'"^) for hard-disk molecules and both of the steady-state BGK equation 
and information theory, respectively. 
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TABLE I: Numerical constants bir in ea. pifl 



rr 


< 1 


r < 2 


r < 3 


r < 4 


r < 5 




r < 6 




r < 7 




1 


1 


1.026 


1.029 


1.030 


1.030 




1.030 




1.030 




2 


- 


5.263 X 10" 


-2 5.657 X 10-2 


5.720 X 10-2 


5.734 X 10- 


-2 


5.738 X 10- 


-2 


5.739 X 10" 


-2 


3 


- 


- 


4.335 X 10-3 


4.820 X 10-3 


4.920 X 10" 


-3 


4.946 X 10" 


-3 


4.954 X 10" 


-3 


4 








3.671 X 10-4 


4.188 X 10^ 


-4 


4.313 X 10" 


-4 


4.349 X 10" 


-4 


5 










2.966 X 10" 


-5 


3.452 X 10" 


-5 


3.583 X 10" 


-5 


6 














2.241 X 10" 


-6 


2.651 X 10" 


-6 


7 


















1.576 X 10" 


-7 
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TABLE II: Numerical constants 6or in ea. H36() 



r r < 2 


r < 3 


r < 4 


r < 5 


r < 6 




2 2.778 


2.530 


2.497 


2.490 


2.488 




3 - 


-3.292 X IQ-^ 


-3.602 X 10-^ 


-3.663 X 10-1 


-3.679 X 10- 


-1 


4 - 




-2.676 X 10-2 


-3.038 X 10-2 


-3.125 X 10" 


-2 


5 - 






-2.294 X 10-3 


-2.670 X 10- 


-3 


6 - 








-1.882 X 10" 


-4 



45 



TABLE III: Numerical constants (upper) and 6^ (lower) in ea. ()37|) 



r 


r < 2 


r < 3 


r < 4 




r < 5 




r < 6 







-1.277 X 10-2 


-1.784 X 10-3 


2.073 X 10- 




6.982 X 10- 


-4 


8.380 X 10- 


-4 


1 


-1.007 


-9.775 X 10-1 


-9.721 X 10 


-1 


-9.708 X 10 


-1 


-9.704 X 10 


-1 


2 


3.363 X 10-1 


3.701 X 10-1 


3.758 X 10- 




3.772 X 10- 


■1 


3.775 X 10- 


-1 


3 




3.262 X 10-2 


3.643 X 10- 


'2 


3.729 X 10- 


2 


3.753 X 10- 


-2 


4 






2.584 X 10- 




2.986 X 10- 


3 


3.089 X 10- 


-3 


5 










2.091 X 10- 


-4 


2.461 X 10- 


-4 


6 














1.557 X 10- 


-5 



r r < 2 r<3 r<4 r<5 r<6 






4.052 X 10-1 


4.041 X 10-1 


4.039 X 10- 


1 


4.038 X 10- 


1 


4.038 X 10- 


1 


1 


-4.261 X 10-1 


-4.297 X 10-1 


-4.304 X 10 




-4.306 X 10- 


-1 


-4.306 X 10- 


-1 


2 


-4.452 X 10-2 


-4.886 X 10-2 


-4.968 X 10 


-2 


-4.987 X 10" 


-2 


-4.993 X 10" 


-2 


3 




-4.594 X 10-3 


-5.174 X 10 


-3 


-5.304 X 10" 


-3 


-5.340 X 10" 


-3 


4 






-4.315 X 10- 


-4 


-4.961 X 10- 


-4 


-5.125 X 10- 


-4 


5 










-3.688 X 10 


-5 


-4.311 X 10 


-5 


6 














-2.881 X 10- 


-6 
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TABLE IV: Numerical constants b2r in ea. H39() 



r 


r < 2 


r < 3 


r < 4 


r < 5 


r < 6 







-2.154 X 10-1 


-2.038 X 10-1 


-2.017 X 10-1 


-2.012 X 10-1 


-2.011 X 10 


-1 


1 


-7.943 X 10-1 


-7.626 X 10-1 


-7.569 X 10-1 


-7.555 X 10-1 


-7.551 X 10 


-1 


2 


3.586 X 10-1 


3.945 X 10-1 


4.006 X 10-1 


4.021 X 10-1 


4.025 X 10" 


■1 


3 




3.492 X 10-2 


3.902 X 10-2 


3.994 X 10-2 


4.020 X 10- 


■2 


4 






2.800 X 10-3 


3.234 X 10-3 


3.346 X 10- 


3 


5 








2.276 X 10-^ 


2.677 X 10- 


-4 


6 










1.701 X 10- 


5 
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TABLE V: Numerical constants for the macroscopic quantities: the ith approximation quantities 
for hard-disk molecules and the exact values for the steady-state BGK equation and information 
theory. 



ith 


Xf 






As 






3th 


-5.085 X 


10" 


-2 


-2.549 X 


10- 


-1 


4th 


-4.807 X 


10" 


-2 


-2.551 X 


10" 


-1 


5th 


-4.757 X 


10- 


-2 


-2.551 X 


10- 


-1 


6th 


-4.744 X 


10" 


-2 


-2.552 X 


10- 


-1 


7th 


-4.741 X 


10- 


-2 


-2.552 X 


10- 


-1 


BGK equation 









1 

4 






information 


1 

2 






1 
4 
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TABLE VI: Numerical constants in ea. ()68|) . 



r 


Boltzmann Eq. BGK Eq. and Information Theory 





-2.758 X 10-2 


3 
128 


1 


-1.761 X IQ-^ 


9 

64 


2 


4.015 X 10"^ 


9 

32 


3 


-1.310 X IQ-^ 


1 

16 


4 


1.326 X 10^2 


- 


5 


-1.368 X 10"^ 


- 


6 


1.260 X lO"'' 


- 


7 


-9.809 X IQ-^ 


- 


8 


6.290 X lO"'^ 


— 


9 


-3.249 X 10"® 


- 


10 


1.306 X 10"^ 


- 


11 


-3.902 X 10-^^ 




12 


8.256 X 10-13 




13 


-1.160 X 10"^^ 




14 


9.645 X 10-^'^ 




15 


-3.576 X 10-19 
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TABLE VII: Numerical constants f3r in ea. (|B^ . 



r Boltzmann Eq. BGK Eq. Information Theory 






-1.284 X 10-^ 


9 

64 


17 
128 


1 


-4.721 X 10-^ 


30 
64 


31 
64 


2 


3.372 X 10-1 


3 
16 


11 
32 


3 


7.086 X 10-2 


1 

8 


1 

16 


4 


-3.874 X 10-3 






5 


1.536 X 10-^ 






6 


-2.774 X 10-*^ 
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